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Abstract 

The question of the value of the pitch angle of DNA is visited from the perspective of a geomet- 
rical analysis of transcription. It is suggested that for transcription to be possible, the pitch angle 
of B-DNA must be smaller than the angle of zero-twist. At the zero-twist angle the double helix is 
maximally rotated and its strain-twist coupling vanishes. A numerical estimate of the pitch angle for 
B-DNA based on differential geometry is compared with numbers obtained from existing empirical 
data. The crystallographic studies shows that the pitch angle is approximately 38° , less than the cor- 
responding zero-twist angle of 41.8°, which is consistent with the suggested principle for transcription. 



In this paper, the restrictions inferred on the double-stranded DNA are studied in a geometrical 
analysis of the mechanics of transcription. Transcription of DNA is the biological process where the 
base-pair code within a gene is read for the purpose of protein expression. The precise mechanisms that 
control transcription are still to be fully explored, and represent a fundamental theme of research in 
molecular biology. In eucaryotes, one fundamental control point for the regulation of transcription is 
the structural adaptation of the chromatin fiber necessary to access the gene sequence. One scenario 
that has been proposed is based on the reversome, explaining how transcription elongation can proceed 
within condensed chromatin 1 . The topology of DNA in the chromosome has been further studied in 
models taking into account possible local modifications and global constraints [2j [3]. For reviews on 
transcription and transcriptional regulation in relation to epigenetic phenomena, see [H [5] . It has early 
been understood that the topology of DNA is essential to understanding transcription, in prokaryotes, 
one suggestion for this is the possible supercoiling of DNA during transcription [HI [7] . 

The mechanics of DNA and its topological supercoiling has been described using elastic models [HI 
[S]. An elastic model has also been used to model DNA loops which play a role for the mechanics of 
transcription |10| ; one suggestion is that transcription of a certain gene can be repressed or activated 
whenever a DNA loop is formed that contains that gene [TT]. In the following, we use an entirely 
geometrical approach and model DNA as a double helix of two flexible tubes with fixed thickness. Under 
the presented assumptions, our main result is that for transcription to be possible, the pitch angle of the 
double helix must be constrained to a certain range, based on a differential length argument. Constraints 
that involve the degree of twisting can have two origins. One origin is of topological nature, such as 
conservation of the linking number [T^] . Another origin is longer-range interactions that are due to the 
range of the involved forces in biological helices jT3l [14] . 

An aspect of transcription involving the pitch angle is the various stresses related to twist and strain 
when the B-DNA reorganizes itself prior to transcription. Experimental studies using magnetic tweezers 
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and optical trapping suggest that DNA has a negative strain- twist coupling [T3 [T^] . Elongated DNA 
will therefore tend to rotate through a larger angle per set of base-pairs, i.e. to be winding-up. Recently, 
an analysis of the observed stick-slip melting of DNA under tension has revealed a zig-zag-like dynamics 
[TT] . In Ref. [T5] , we have suggested a geometrical model that describes the phenomenon of winding-up 
of biological double helices as a generic phenomenon for all molecular double helices for which the pitch 
angle is below a specific value. 

In the following, we review the geometry used in the analysis and study its implication on the pitch 
angle. Later, we compare with the experimentally available data and determine the pitch angle of B-DNA. 
Twisting two strands together to form e.g. a geometry akin to the double-stranded DNA, the resulting 
double helical structure has a length, L r (n), that depends on the number of twists, n [TS]. When the 
strands are parallel to each other, the length of the double helix is the same as the length of the strands, 
and the double helix becomes shorter and shorter as turns are added. Only up to a maximum number 
of turns can be added, resulting ultimately in a geometry which is maximally rotated. In addition, this 
structure has a vanishing strain-twist coupling |18j . As rotations are removed from the helix, the length 
becomes further reduced. Therefore, L r (n) is shaped like a horizontal hairpin, see Fig. [I] 

In the detailed calculation of L r (n), the double helix is simplified to be consisting of two circular 
tubes of diameter D, whose center lines are coiled with pitch H on a common cylinder of radius a. The 
pitch angle, v±, measured from the horizontal is determined by t&nv± = H/2i:a. For these idealized 
structures, the steric interactions of molecular DNA are described by letting the tubes be in contact with 
hard walls, i.e. the distance between the two helical lines should satisfy 



min(|ri(ii)-r 2 (i 2 )|) =D, (1.1) 

where ft, r 2 describe the two helical center lines. 

For transcription we consider DNA to be part of a structure that does not change on a large scale. 
Within this structure we consider an intermediate structure which is allowed to change while being 
maintained as a double helix. Finally, we consider a relatively short stretch of DNA that separates 
(melt) to single-stranded DNA. The first of these considerations require that the length of DNA during 
transcription, L' DNA , to be greater than the length of DNA before transcription, Ldna, i-e. 

L'dna>^dna- (1-2) 
Now, the twisting behavior of DNA depends on its pitch angle, v±. If the pitch angle is such that it 



would just unwind under strain, then Eq. (1.2) would not be fulfilled, as shown below. Instead the DNA 
would tend to get shorter during transcription. This is the origin of the constraint on v± and shows that 
the behavior of B-DNA under strain must be opposite of that of a rope. 



Mathematically, the two helical center lines in Eq. (1.1) have the parametric equations 



^i(^i) = (acosfi, asinti, (H/2ir)ti) , 

r 2 {t 2 ) = (acosi 2 ,asint 2 ,(#/27r)i 2 + *(#/27r)), (1.3) 

with £x , £ 2 € R" The parameter ^ is a phase parameter that determines if the double helix is symmetric 
or asymmetric. For a symmetric double helix, e.g. A-DNA, ^ is equal to 7r and for an asymmetric double 
helix, akin to B-DNA, we have * = 7rl44°/180° = 2.513. The phase difference of 144° for the minor 
groove in B-DNA is described in Ref. [20]. With D being the tube diameter, the points of tube-tube 
contacts obey the equation 
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Figure 1: Relative length, e, of double helix as a function of the relative number of turns, t, on the strand. 
The curve is plotted for the asymmetric case, '5 = 2.513. The solid curve is the configurations where the 
two strands are in contact. The pitch angle v± increases along the curve from 0° at the lower left corner, 
to 41.8° at the turning point, and finally to 90° at the upper left end-point. In the present paper B-DNA 
is found to have a pitch angle of ~ 38° corresponding to the point below the turning point. The contour 
lines have constant volume fraction, along the outer contour its value is 0.597 corresponding to B-DNA. 



H ,2 



D 2 = a 2 {cost- l) 2 + a 2 sin 2 i + ( — ) (* + t) 2 , (1.4) 



where t = t\ — t%. The solutions to the condition of minimum distance, Eq. (1.1 1, is found from a 
numerical study of the equation, dD 2 /dt = 0; for further details, see Refs. [2"Tl 1221 12*3] . These solutions 
have recently been expressed in a non-Euclidian geometry |24j . 

Figure [l] shows how the (relative) length of the double helix depends on the number of turns in the 
twisted strands using Eq. \1A\ ; the specific curve is for the asymmetric double helix with "J = 2.513. 
The vertical axis is the relative length, e, of the double helix, where e = L r /L s and L s is the strand 
length. The horizontal axis is the relative number of turns, r = n/n max , where n max is the number of 
turns at the turning point for L s — 1. At the turning point, where r is maximal, the pitch angle of the 
double helix takes a unique value, v± = vzt, at which the double helix has zero strain- twist coupling. In 
[T8] it was found that vzt = 41.8° for an asymmetric double helix (B-DNA) corresponding to e = 0.6665. 
The contour plot in Fig. [T] shows lines of constant volume fraction, fy. The volume fraction is defined 
locally as the volume of the two strands relative to the volume of a smallest enclosing cylinder [21 . The 
boundary of the contour plot has fy — 0.597 which is the maximal volume fraction for the asymmetric 
double helix. The contour plot shows that when straining DNA, the geometry will tend to follow the 
solid curve as its tangent also is tangent to lines of constant volume fraction. 

A double helix structure with a pitch angle greater than the zero- twist angle, v± > vzt, is at the upper 
part of the curve in Fig.[l] Now, there are two possibilities for the DNA double helix. Either v± > vzt, or 
vj_ < vzt- We will argue that actually, v±_ < vzt, so that DNA cannot be a maximally twisted structure, 
and that this bound is in agreement with the principle of close-packing of helices [21] . Let us assume that 
transcription is to take place on a short stretch of a long stretch of DNA without changing the total twist. 
The long DNA has strand length L, and the short stretch strand length I. The number of turns on the 
long DNA is u^na = tl ■ n max ■ L, and during transcription it is n' DNA = Tl—i ■ n max ■ (L — I). Therefore, 
tl-i = On the stretch where transcription takes place, we have e; = 1; that corresponds to the 

upper left endpoint of the curve in Fig. [l] The remaining part of the DNA has a relative length e^-i, 
which for now is undetermined. The difference in the length of DNA after and before transcription is, 

l 'dna - Ldna =1- ei + (L-l) ■ e L -i - L ■ e L . (1.5) 
If this number is positive, then the DNA gets longer under transcription. If it is negative the DNA gets 
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Figure 2: B-DNA structure PDB 425D 



shorter, which is not generally physically possible. We assume that I is small compared to L, I <C L, and 
series expand cl-i to obtain, e^_; = e; + ^ i~i t l- Therefore Eq. (1.5) can be rewritten as 



L'dna ~ Ldna =£11 



de 



-T~L 



(1.6) 



Due to the convex nature of the curve in Fig. [T] the requirement that this difference is greater than zero 
cannot be obtained on the curves upper part where it is continuously bending downwards. We must be 
on the lower part of the curve in Fig. [I] so that the pitch angle obeys v± < vzt- 

What is the experimentally determined pitch angle, v±, of B-DNA? There is some debate in the 
scientific community about the interpretation and specific value of the pitch angle of B-DNA. In particular, 
some publications advocate a structure for DNA that has a pitch angle of 45° [25j [26] based on the 
symmetric double helix. This corresponds to a point on the upper part of the A-DNA symmetric version 
of the curve in Fig. [T] with e = 0.707. We will argue that the pitch angle for B-DNA is about 38°. We 
determine the pitch angle of DNA in the following way. The DNA pitch height, H, can accurately be 
determined by experiment. It is the height for which the helical structure progresses with precisely one 
full rotation. The helix-packing diameter [27], or simply the width, W, of the molecule, can also be 
accurately determined. From H and W the aspect ratio can be calculated as, 



A = H/W. (1.7) 

A classical work that lists the values of H and W for crystal structures of DNA is the Cold Harbor 
Symposium contribution by Dickerson et al. 27J. The results are reproduced in Table [TJ together with 
our calculated aspect ratios and the obtained pitch angles. Figure [2] depicts H and W for the B-DNA 
structure PDB 425D [35]. The helical geometry is clearly visible, however it is superficial to determine 
the pitch angle by the naked eye from a two-dimensional projection. One needs to carry out a genuine 
three-dimensional geometrical study. 

The pitch angle can be indirectly determined from the aspect ratio above. For the double helix with 



phase parameter 'J', the equation ( 1.4) gives a unique curve for 2a/ D as a function of pitch angle, that is 



monotonically decreasing. The relation between the pitch angle v± and the aspect ratio A reads, 

. . . 7rtanwj_ ,„ „ N 

A(v ± )= , 1.8 

<p{v±) 
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Figure 3: Aspect ratio A versus pitch angle v± plotted for double helix with minor groove (solid, red curve) 
and symmetric double helix (dashed, blue curve). The solid circles are for A-DNA, and B-DNA using 
Ref. JZTf. 



where (f){v±_) = 1 + D/2a. Figure [3] shows how the aspect ratio, A, depends on the pitch angle, vj_, for 
a tubular double helix by solving numerically Eq. (1.1 1. The solid curve in Fig. [3] is for an asymmetric 



double helix mimicking B-DNA. For an aspect ratio of 1.40 we obtain the pitch angle v± = 38° shown by 
a solid circle. The dashed curve is for a symmetric double helix such as A-DNA. It passes through the 
point (45°,7r/2). The experimentally determined aspect ratio for A-DNA is 0.96. This corresponds to a 
pitch angle of 28° as shown in Fig. [3] by a second solid circle. 

In its crystallized form, B-DNA has 10 base-pairs per turn, and using Eq. ( |1.8[ ) it gives a result for 
the pitch angle which is ~ 38°. Non-crystallized DNA, on the other hand, has about 10.5 base-pairs per 
turn. It is therefore further away from the turning point of the curve in Fig. [T] and consequently has a 
slightly lower pitch angle [53]. Table 1 summarizes our results for crystallized DNA. 
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B-DNA 


33.2 


23.7 


1.40 


38° 41.8° 


A-DNA 


24.6 


25.5 


0.96 


28° 39.4° 



Table 1: Pitch height, H , and width, W , experimental values for crystallized A-DNA and B-DNA re- 
produced from Ref. \2T$ . Also given are derived values for aspect ratio A and pitch angle v±, and the 
corresponding zero-twist pitch angle vzt from Ref. U8f . 

The obtained pitch angle for DNA is therefore in agreement with the bound determined from transcription. 
It can however seem a surprising result, how can the DNA pitch angle be below 45°? In the symmetric 
double helix of ^ = tt, the two strands have a straight line of contact for v± > 45° [30 . Outside this 
interval, the two strands have a helical line of contact and the double helix has a central channel. In 
Ref. [53] , an elastic model is discussed where the pitch angle is to be greater than or equal to 45° , which is 
suggested to be a lock-up angle for a two-ply [3] ■ For the purpose of predicting the molecular interactions 
underlying the value of the pitch angle of DNA an elastic model with energy terms for curvature and 
torsion is not in agreement with the proposed bound. Presumably such an elastic model does not describe 
well the interactions at play on the molecular level. 



Both of the empirical numbers for v± are in good agreement with the criterion suggested in |21] that 



the dominating packing feature of DNA is determined by the obtainable volume fraction. I.e. that the 
criterion is consistent with an optimum use of available space. The packing of helical DNA therefore 
follows a classic packing rational in condensed matter [31] . It seems pertinent to pursue optimal packing 
in terms of maximization of the volume fraction, and hence density. The theoretical optimal values using 
the volume-packing criterion is 32.5° for the symmetric helix J3TJ[35] an d 38.3° for the asymmetric helix 
[2"T] . both less than the corresponding zero-twist angles so that transcription can take place. 
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